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Abstract 

Sources of gravitational waves (GW) and emitters of high energy (HE) neutrinos both involve compact objects and 
matter moving at relativistic speeds. GW emission requires a departure from spherical symmetry, which is the case 
if clumps of matter are accreted around black holes or neutron stars, and ejected in relativistic jets, where neutrinos 
are believed to be produced. Both messengers interact weakly with the surrounding matter, hence point directly to 
the heart of the engines that power these emissions. Coincidences between GW interferometers (e.g. Virgo) and HE 
v telescopes {e.g. Antares) would then give a unique insight on the physics of the most powerful objects in the 
Universe. The possibility, observability and detectability for such GW/HE v coincidences are analysed. 
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1. Scientific motivations 

The forthcoming years should be very excit- 
ing both in gravitational wave (GW) astronomy 
and high energy neutrino (HE v) astronomy. The 
ViRGO interferometer [1] , currently closed down for 
upgrade, should be taking data with an improved 
injection system in 2009, whereas the Antares col- 
laboration has completed the deployment and 
connections of its 12 lines [2] early in June 2008, 
starting the operation of the first undersea neutrino 
telescope. In the mean time, the two LiGO interfer- 
ometers (ITF) are in operation [3J, and IceCube 
is deploying its 1 km 3 neutrino telescope in the ice 
of the South Pole, already having 40 lines in data 
taking mode [4] . 

Both GW sources and HE v emitters involve 
compact objects and matter moving at relativistic 
speeds. As a result, coincidences between GW in- 



terferometers and neutrino telescopes can be envis- 
aged [5J, and should then be feasible by 2009, with 
the joint operation of Virgo, Ligo, Antares and 
IceCube. Together, they would give a unique in- 
sight on the physics of the most powerful objects 
in the Universe. Some classes of astrophysical ob- 
ject, invisible in electromagnetic channels, may be 
observable only via their gravitational waves and 
high energy neutrino emissions. Finally, in many 
quantum gravity (QG) models ,6], the propagation 
velocity of a particle depends on the energy: 
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where |£| ~ 1, and Eq G is the energy scale at which 
QG effects arise. Hence, measuring a non-zero time 
delay between gravitational wave bursts and high 
energy neutrino transients would then allow to probe 
Quantum Gravity effects at the Planck energy level. 
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Such time coincidences require GW bursts (local- 
ized in time) , and would allow GW antennae to con- 
firm a burst detection, and neutrino telescopes to 
sign the cosmic origin of the signal. It must be first 
demonstrated that such a coincidence process is: 

• possible: sources of both GW and HE v, able to 
emit signals possibly coincident in time, exist in 
the Galaxy. This is discussed in section [2] 

• observable: the visibility sky maps of ViRGO and 
Antares are not orthogonal. This point is ad- 
dressed in section [3] ; 

• detectable: the resulting coincidence detection 
probability, for a fixed accidental coincidence 
rate, which constrains the authorized background 
level for each individual detector, must be esti- 
mated. This point is developed in section [4] 

2. GW bursters and HE v sources 



time scale of the event. The duration of major ejec- 
tion events is typically on the order of days, while 
that of less powerful flares is correspondingly shorter 
(minutes to hours). The correlations between the 
X-ray and synchrotron emission clearly indicates a 
connection between the accretion process and the jet 
activity. Whether radio and IR outbursts represent 
actual ejection of blobs of plasma or, alternatively, 
formation of internal shocks in a quasi-steady jet is 
unclear. In any case, since the overall time scale of 
outbursts is much longer than the dynamical time 
of the compact object (milliseconds), it is likely that 
shocks will continuously form during the ejection 
event. If a fraction of at least a few percent of the 
jet power is used to accelerate electrons to very high 
energies then emission of high-energy gamma rays 
is anticipated, in addition to the observed radio and 
IR emission. 



Only galactic potential sources of both GW and 
HE v signals will be discussed here, for two reasons. 
Firstly, the only sources accessible with first gener- 
ations detectors such as Antares and Virgo are 
likely to be galactic ones. Secondly, for sources in 
the local universe, in particular with negligible red- 
shift, the QG delay mentionned in the previous sec- 
tion is independent of cosmological models. This 
(non-exhaustive) section focuses on microquasars 
and soft-gamma repeaters, however core-collapse su- 
pernovae and gamma-ray bursts are commonly cited 
extra-galactic sources of both GW and v emissions: 
these are described elsewhere [7]. 

2.1. Outbursts from microquasars 

Microquasars are galactic jet sources associated 
with some classes of X-ray binaries involving both 
neutron stars and black hole candidates. During ac- 
tive states, the X-ray flux and spectrum can vary 
substantially, with a total luminosity that, during 
the so-called very high states, often exceeds the Ed- 
dington limit [5]. Their activity involves ejection 
within jets with kinetic power that appears to con- 
stitute a considerable fraction of the liberated ac- 
cretion energy, giving rise to intense radio and IR 
flares. Radio monitoring of some X-ray transients 
has revealed supcrluminal motions in some objects, 
indicating that the jets are relativistic, with 7^1 — 
10. The ejection episodes are classified into several 
classes according to the brightness of synchrotron 
emission produced in the jet and the characteristic 



2.1.1. High energy neutrino emission 

The content of jets in microquasars remains an 
open issue. In scenarios in which an initial rise of the 
X-ray flux leads an ejection of the inner part of the 
accretion disk, as is widely claimed to be suggested 
by the anticorrelation between the X-ray and radio 
flares seen during major ejection events, e-p jets are 
expected to be produced. A possible diagnostic of 
e-p jets is the presence of Doppler-shifted spectral 
lines, such as the Ha line as seen in SS433. 



z(production) = 10 cm in jet 
{base of jet) 
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Fig. 1. 7 Data from EGRET and HESS, together with pre- 
dictions for models with production of 7 at the base or at 
the end of the jet. 

Taking the example of LS 5039, some authors ar- 
gue in favor of hadronic origin of TeV photons |3] , 
especially if they are produced within the binary 



system. The detected 7-rays should be accompanied 
by a flux of high energy neutrinos emerging from 
the decays of ^ mesons produced at pp and/or 
P7 interactions. The flux of TeV neutrinos, which 
can be estimated on the basis of the detected TeV 
7-ray flux, taking into account the internal 77 — > 
e + e~ absorption, depends significantly on the loca- 
tion of 7-ray production region. The minimum neu- 
trino flux above 1 TeV is expected to be at the level 
of 10 12 cmV 1 , but could be up to a factor of 100 
higher. As can be seen in figure[l](taken from [9]) , the 
HESS/EGRET data agree well with a production of 
7 (and neutrinos) at the base of the jet, very close 
to the onset of the acceleration phase and its corre- 



sponding GW signal (see section 2.1.21. Finally, one 



should note that the detectability by Antares or 
future km 3 telescopes depends strongly on the high 
energy cutoff in the spectrum of parent protons. 
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2.1.2. GW emission during a powerful flaring event 
Two kinds of processes could lead to detectable 
GW signals [10]. Firstly, the matter accreted for 
months/years could be swallowed by the compact 
object, and, provided that the process is fast, trigger 
the resonance of normal modes in the central object. 
Secondly, the acceleration of the matter in the jet is 
the origin of a short GW burst. 

For both signals, the amplitude will depend criti- 
cally on the accreted/ejected mass. The assumption 
which is made in the following is that all the matter 
ejected during a flare comes from the accretion disk, 
and had previously fallen onto the compact object 
at some previous moment. The GW emission pro- 
duced by the infall of matter onto the central ob- 
ject is different in the case of neutron stars (NS) and 
black holes (BH) , but in both cases results in the ex- 
citation of the Quasi-Normal Modes of the star [11] 
(typically a damped sine signal) which could con- 
tinue into the ejection phase. The time-lag between 
the two processes is unknown, and could range from 
ms up to several days. 

In the NS case, the characteristic amplitude 
(scaling as 1/d) at 1 kpc can be written as 

1 icy 1 

h c « 10~ 20 (xjj^e) /o 2 for the accretion sig- 
nal [11] , to be compared with the ViRGO sensitivity 
/inoisc ~ 10~ 21 around 1 kHz. 5m — r/M© is the 
accreted/ejected mass, and / the frequency of the 
excitation. Note that this assumes a sudden infall 
of the matter onto the compact object (bursting 
mode) , and does not stand if the accretion infall is 
a continuous and slow process (continuous mode). 



Fig. 2. Amplitudes expected for the accretion/ejection pro- 
cesses at 1 kpc, for 5m ~ 1O _4 M0, race ~ 10 ms, 7 = 10. 



The acceleration of an ultrarelativistic blob of 
matter with a Lorentz factor 7 around a compact 
object induces a burst with memory, the space-time 
perturbation reaching a maximum amplitude at 1 
kpc Ah ~ 2x 10~ 22 (jq) (j^e), independent on 
the density of the ejecta [12]. The frequency is typi- 
cally the inverse of the acceleration time £ acc , rang- 
ing from A< fl . oc _f all ~ 0.1 ms for a 10 M© object, up 
to Ar max ~ 1 min, observations of galactic X-Ray 
binaries showing that radio emission occurs at a dis- 
tance ~ 0.1 AU from the central object. A summary 
of those estimates are displayed in figure [2] for a 
source at 1 kpc, ejected mass 5m ~ 10~ 4 M and 
acceleration time r acc ~ 10 ms. 

Taking into account accretion rates (typically 
M ~ 10~ 8 Af Q /yr), jet luminosities (10 29 J/s - 
10 33 J/s) and radio blob sizes ~ 10 12 m and den- 
sities (from 10~ 10 g.cm~ 3 up to 10 15 g.cm -3 ), the 
accreted/ejected mass for major events can be esti- 
mated to range from 10 _8 M Q up to 1O~ 4 M0, with 
an average mass of 10 _6 M Q . Thus the amplitudes 
shown in figure [2] correspond to the maximum mass 
and Lorentz factor that can be expected. Nonethe- 
less, such extreme outbursts, where the matter ac- 
creted over months or even years, is swallowed and 
then ejected in one single bursting event cannot be 
ruled out. 



2.2. Flares from Soft-Gamma Repeaters 



3. Observability of Coincidences 



Soft gamma-ray repeaters (SGRs) are X-ray pul- 
sars which have quiescent soft (2-10 keV) periodic 
X-ray emissions with periods ranging from 5 to 10 s. 
They exhibit repetitive bursts lasting ~ 0.1 s which 
reach peak luminosities of ~ 10 34 J/s, in X-rays and 
7-rays. There are 4 known SGRs, 3 in the Milky Way, 
and one in the Large Magellanic Cloud. Three of the 
4 have had hard spectrum (MeV energy) giant flares 
with luminosities up to 10 40 J/s. The favoured mag- 
netar model for these objects is a neutron star with 
a huge magnetic flied B ~ 10 15 G [13]. Star-quakes 
are thought to fracture the rigid crust causing out- 
bursts. These giant flares result from the formation 
and dissipation of strong localized currents coming 
from magnetic field rearrangements associated with 
the quakes, liberating a high flux of X- and 7-rays. 

Sudden changes in the large magnetic fields would 
accelerate protons or nuclei that produce neutral 
and charged pions in interactions with thermal ra- 
diation. These subsequently decay into TeV or even 
PeV energies 7-rays and neutrinos |14j . Flares from 
SGRs are thus potential sources of high energy neu- 
trinos. 

During the crustal disruption, a fraction of the 
initial magnetic energy is annihilated and released 
as photons, and the stored elastic energy is also con- 
verted into shear vibrations with frequencies in the 
kHz regime. These waves are able to excite non- 
radial modes, damped by GW emission [T5]. The ex- 
pected gravitational strain amplitude can be writ- 
ten as: 
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where E is the total energy, G the gravitational con- 
stant, Lo n and r„ the pulsation and damping time- 
sale of the n-mode. These oscillation parameters de- 
pends on the equation of state and the stellar mass. 
Simulations |15j show that low mass stars produce 
larger gravitational amplitudes and could be de- 
tected more deeply within the Galaxy. Moreover, 
sources at distances ranging from 0.4 kpc up to 2.4 
kpc could be probed with the planned sensitivity of 
ViRGO. The detection probability and frequency of 
these events depend on the poorly constrained dis- 
tribution of this class of sources in the Galaxy. 



Both GW detectors and HE v telescopes have 
limited sky coverage and exposure. In order to per- 
form coincidences between both types of detector, 
the overlap of such visibility maps has to be com- 
puted: the following paragraphs address this ques- 
tion, taking the examples of ViRGO and Antares. 

3.1. ViRGO beam pattern 

The response h(t) of an interferometric detector 
to a GW is a linear combination of the two indepen- 
dent wave polarizations h + and h x , with weighting 
factors called the beam pattern functions: they have 
values in the range [—1; 1], depending on the longi- 
tude and latitude of the detector location, as well as 
its orientation, the angle between the arms, the sky 
coordinates of the source, and the wave polarization 
angle. The best response is achieved for detectors 
with orthogonal arms. In the following, average over 
the unknown polarization angle will be presented. 
The instantaneous beam pattern (normalised to its 
maximal value) at a given time during the day, in 
equatorial coordinates (right ascension a, declina- 
tion S), is displayed in figure |3]for ViRGO. 
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Fig. 3. ViRGO beam pattern in equatorial coordinates. 



3.2. Antares visible sky 

Antares is only sensitive to sources below the 
horizon at some time during the day, because it 
searches for neutrinos which have interacted in the 
Earth: a portion of the sky is therefore never visi- 
ble. Figure [4] shows the daily average visibility as a 
function of sin 8. 
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Fig. 4. Antares daily averaged visibility vs sin 8. 

3.3. Antares/Virgo common sky 

The visibility sky map for coincidences between 
Antares and Virgo is the convolution of the two 
previous exposure maps. The daily averaged com- 
mon sky map is displayed in figure [5] together with 
the position of known microquasars and soft-gamma 
repeaters (or magnetars). Except for three of them, 
all are visible at some time by both experiments, ren- 
dering observable GW/HE v coincidences for most 
of these galactic sources. 

■ MicroQuasars ^Magnetars (SGRs) 
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Fig. 5. Common visibility sky map for Antares and Virgo. 



4. Detectability of Coincidences 

Given GW/HE v coincidences are observable, at 
least for the Antares and Virgo couple, the de- 
tectability, in terms of detection probability vs false 
alarms rate, has to be evaluated. 

4.1. Coincidence window and QG effects 

To set the coincidence time window, possible 
physical propagation delays have to be estimated. 



In the case of GW, the graviton being massless, and 
the energy carried away by each individual graviton 
in a GW burst being small (-E grav iton ~V<1 for 
/ = 1 kHz), QG-induced or mass-induced delays 
are close to zero. For a 1 TeV v, the mass-induced 
delay is negligible even with m„ = 1 cV. Taking the 
expression given in eq. [T] as a starting point, with 
£ = — 1 being favoured, and neglecting any cosmo- 
logical effects (for low redshift z <C 1), the delay in 
ms becomes, in the first order: 
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Taking E QG — -E P i anck ~ 10 19 GeV, this yields a 
maximum QG delay of 1 second for E v = 1 PcV 
and sources up to the Large Magellanic Cloud (d ~ 
50 kpc) , or lower energy neutrinos (E v w 1 TeV) and 
sources as far as the Virgo Cluster (d ~ 20 Mpc). 
At coinc = 1 s thus seems a reasonnable choice. Never- 
theless, the coincidence time window can also be set 
by imposing an overall coincidence detection proba- 
bility for a given GW signal. Detection issues, both 
for GW ITFs and HE v telescopes must now be ad- 
dressed. 



4.2. Virgo detection 

The detection of a transient signal in a GW ITF is 
not an easy task. Generally, an allowed false alarm 
rate is fixed, and the detection probability can be 
estimated as a function of the signal-to-noise ratio 
(SNR or p) of a particular signal: this is shown for 
ViRGO in figure [6] for p max = 5 (defined as the SNR 
obtained with optimal detector orientation and per- 
fect detection by Wiener filtering), a very low sig- 
nal, for detection algorithms designed for burst sig- 
nals [16 . In the best case, a threshold corresponding 
to 1 false alarm every 5 minutes is needed to obtain 
a 50% detection probability, without taking into ac- 
count any beam pattern effects. 

The right part of figure [6] displays the detection 
efficiency as a function of /9 max , for this particular 
false alarm rate, in the case of a single ITF detec- 
tion or coincident detection in the ViRGO/LiGO net- 
work: for low SNRs, the detection by a single ITF 
is more probable than any twofold coincidence, and 
the detection by any single ITF is always more effi- 
cient than any coincidence configuration |17j . In the 
case of the detection by any of the 3 ITFs, the di- 
rectional information is not available (no triangula- 
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Fig. 6. Left : detection probability vs false alarm rate, for a 
p = 5 gaussian pulse. Right : detection probability vs p ma x 
for the different possible detection configurations. 



tion) , and the only relevant information is therefore 
the time of the burst event. 

ViRGO has a sampling frequency of 20 kHZ, and 
for a gaussian burst of width r and SNR p, the rms 
error on the burst arrival time is |17j : 



Ae 



1.5 

SNR 



VI ms/ 



(4) 



yielding a timing resolution below the ms for p > 5 
and short burst. This of course limits the accessible 
QG energy scale, and the coincidence window to be 
used. 



4.3. Antares detection 

In a neutrino Telescope, the Cerenkov light emit- 
ted by the neutrino-induced muon is detected by an 
array of photomultipliers arranged in strings, able to 
reconstruct the energy and direction of the incident 
muon/neutrino. The measurements of the time of 
the hits and the amplitude of the hits, together with 
the position of the hits are needed to achieve the 
reconstruction of the muon track with the desired 
resolution (below 0.3° above 10 TeV). The quality 
of the track fit is often expressed in terms of a log- 



likelihood ratio term A 



iog(£) 



the distribution 



of which is shown in figure u\ for (upward) atmo- 
spheric neutrinos and misreconstructed atmospheric 
muons, together with the signal detection efficiency 
(i.e. the ability to detect within 1° of the true direc- 
tion a signal neutrino, assuming a E~ 2 spectrum). 
The standard cut applied is A = —5.3, for which the 
signal efficiency is close to 75%; the misreconstruted 
atmospheric muons are decreased to 1/day, and the 
atmospheric neutrinos to 10/day [15] . 



>. 
















?10 4 

1) 

Q. 

£ 3 
C 10 

> 
<u 






an DacKgrouna 
[ upw. rec. atm. it 
atm v 

signal efficiency 










?10 2 

"E 
'(5 

£10 






r | 


: -4 : , "IP 




50% 


I 4 evts/day=p.iu /s 1 








o 
S 1 








%L 










i >. 


E 

2io" 1 

10* 










^ 


=, , , , I , , , , I , , , , i , i , , I , , , , 




, , i 


..l..M 



-8 -7.5 



-6.5 -6 



-5.5 -5 -4.5 -4 
value of cut on A 



Fig. 7. Number of background events left as a function of 
the cut value for A. The efficiency for signal is also shown. 

4.4. Accidentals and efficiency 

Figures[6]and|7]provide the information needed to 
estimate the detection probability ey,A for a back- 
ground/false alarm level Rv,a m both detectors (V 
for ViRGO, A for Antares). The coincidence de- 
tection probability is e coinc = eye^, whereas the co- 
incident accidentals rate in a given time window is 
■Rcoinc = RvRAAt coinc . Setting At coinc = 1 s and 
-Rcoinc ~ l/y r i the resulting coincidence detection 
probability is shown in figure [8] as a function of the 
A cut: the efficiency is maximum for A ~ —5.5, be- 
low which Ra is too high, resulting in a too high 
ViRGO detection threshold (for a preset RvRa)- 

Equivalently, the coincidence detection probabil- 
ity can be set at e.g. 50% for a given signal, and 
the maximal allowed coincidence time window can 
be extracted: this is displayed in figure [9] as a func- 
tion of Antares/Virgo detection probabilities, for 
a p = 5 gaussian burst, with i? CO i nc ~ 1/yr. The 
time coincidence window is maximal for ty ~ 65%, 
reaching ~ 15 ms. 

For such a low accidental rate (1/yr), several coin- 
cident detections would be needed to have a high sig- 
nificance, but this nonetheless proves the detectabil- 
ity of GW/HE v coincidences. 

5. Antares/Virgo Coincidences 

If a coincidence is indeed observed, the signifi- 
cance of this positive detection can be estimated, 
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Fig. 8. Coincidence detection probability vs A cut. 
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Fig. 9. Maximal allowed coincidence time window in 
ms vs Antar.es/Virgo detection probabilities, requiring 
e coinc = 50%. 

and assuming that both the GW and HE v signals 
have been emitted with zero delay at the source, 
limits can be put on the QG energy scale -Eqg- I n 
the case described in the previous section (requiring 
£coinc — 50%), the minimum accessible -Eqg can be 
estimated using eq. [T] The maximum energy scale 
yielding a measurable effect is limited by GW tim- 
ing resolution, which depends on the burst dura- 
tion and SNR, and reaches in this case £"q G x ~ 
5 x 10 18 GeV (p = 5), close to the Planck limit. The 
minimum accessible energy scale is in turn deter- 
mined by the maximal coincidence window defined 
previously, which yields Eq£ ~ fO 17 GeV. This is 
to be compared with existing limits on Eqq, e.g. us- 
ing TeV flares from Mrk42f ~ 4 x fO 17 GeV [19]. 



It should be noted that to perform a real measure- 
ment of Eq G , the neutrino energy resolution is of 
importance, and is a factor 2 or 3 in the case of 
Antares H"8l. 



5.1. Making coincidences 

The process for performing time coincidences 
is quite classical. Virgo and Antares produce 
trigger lists, according to some predefined false 
alarm rates, corresponding to different coincidence 
time windows. Timeshifts performed on these data 
streams allow for the study of background coincident 
triggers; in the zero-lag case, Antares/ ViRGO trig- 
gers are compared in predefined time windows. 
The significance of an observed coincidence is then 
statistically evaluated by comparing the two cases. 



5.2. Antares/Virgo common calendar 

ViRGO took data jointly with the 2 LlGO interfer- 
ometers between May and September 2007, during 
the Virgo Scientific Run (VSR), achieving the sen- 
sitivity shown in figure |10[ together with the sensi- 
tivities obtained during previous commissioning or 
science runs (weekly runs from September 2006 un- 
til march 2007) . In spite of problems with the laser 
injection system (at low frequency), and a factor 2 
difference at high frequency, the resulting sensitiv- 
ity is quite close to the expectation. The interferom- 
eter should be taking data again with an improved 
injection system in 2009. 

Antares has been taking data continuously with 
its final 12 line configuration since the end of May 
2008, and is expected to observe high energy neu- 
trinos for a period of 10 years. Interestingly, during 
the VSR, Antares already had 5 lines operational 
since January 2007. Clearly, during this parallel op- 
eration of Virgo and Antares 5-lines, only the 
most powerful GW/HE v sources could be detected; 
however these data could be used as a test bench for 
preliminary studies on time coincidences. Moreover, 
LlGO data, which are also available for this VSR/5- 
lines period, could be used to enhance the detection 
probability, a detection by any of the 3 interferome- 
ters being more efficient than for ViRGO alone (see 
section [I| . 
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Fig. 10. VlRGO sensitivity curves from Commissioning Runs 
to Science Runs. 



6. Conclusions 

Time coincidences between GW interferometers, 
such as VlRGO, currently down for upgrade, and 
high energy neutrino telescopes, such as Antares, 
now fully operational, are thus feasible: common 
sources able to emit coincident signals in both chan- 
nels exist in our own galaxy. Microquasars, dur- 
ing major outbursts, and flares from SGRs are pos- 
sible targets. Antares and Virgo visibility sky 
maps are not othogonal, allowing for coincident de- 
tections. Taking into account the relationship be- 
tween detection efficiencies and false alarm rates in 
Antares and Virgo, coincident observations ap- 
pear possible. 

Such coincidences can be performed using 
Antares 5-lines and VSR data (end of 2007) as a 
test that could be performed before the upgrade of 
VlRGO next year, which should improve the sensi- 
tivity at low and high frequency to reach the design 
sensitivity. This VIRGO+ upgrade will correspond 
with the routine operation of the full Antares de- 
tector, provinding the opportunity to perform the 
time coincidences presented in this paper. 

Finally, circa 2015, a km 3 neutrino telescope 
should be operating in the Mediterranean Sea [2U] , 
along with an Advanced VlRGO interferome- 
ter |21j . with enhanced sensitity at low frequency. 
Figure 11 [TO] shows the ejected mass needed in a 
microquasar ejection event to obtain a SNR = 5 in 



Virgo+ and Advanced Virgo [21] as a function 
of the acceleration time (see section[2]): less extremes 
ejection scenarios could be probed, and interesting 
constraints on accretion/ejection models could be 
set by this novel multimessenger approach. 
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Fig. 11. Ejected mass needed for SNR=5 for a microquasar 
at 1 kpc for VIRGO+ and Advanced Virgo. 
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